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Abstract 

Based on the Schwinger- Dyson (SD) equation, the fermion mass generation is 
further studied in the D(2 < D < 4)-dimensional Thirring model as a gauge theory 
previously proposed. By using a certain approximation to the kernel, we analytically 
obtained explicit form of the dynamical mass of fermion and the critical line in (N, 1/g) 
space, where N is the number of fermions and g is the dimensionless vector-type four- 
fermion coupling constant. This analytical result is confirmed by the numerical solution 
for the SD equation with exact form of the kernel in (2+1) dimensions. 
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§1. Introduction 

Through the studies of the dynamical origin of the quarks/leptons mass, it was dis- 
coveredS*'!)'!) that the dynamics inducing the large anomalous dimension of the operator 
ijjijj is strongly required. In particular, the scalar/pseudoscalar-type four-fermion inter- 
actions with the gauge interaction (gauged Nambu-Jona-Lasinio modeltP' have played 
a very important role as a renormalizable model in (3+1) dimensions (D = 4).0 ) It has 
also been shown B that the phase structure of such a gauged NJL model in (3+1) di- 
mensions is quite similar to that of the D{2 < D < 4)-dimensional four-fermion theory of 
scalar/pseudoscalar-type (without gauge interactions), & often called Gross-Neveu mo deli 
which is renormalizable in 1/N expansion. 0) 

Recently a model with only the vector-type four-fermion interaction has been studied as 
another candidate for the scenario of the fermion dynamical mass generation, 0)>El)>lll)>0)> © 
namely, the D(2 < D < 4)-dimensional Thirring model© with N fermions. If we naively 
treated this model with the usual gap equation made by 1/N- leading diagram, we would 
conclude no fermion mass generation. Nevertheless, we can obtain the dynamical mass 
of fermion by solving the SD equation with the rainbow diagram of the composite vector 
boson which is induced in 1/N- leading order. 0* Then the dynamical mass has a nontrivial 
dependence on 1/N as a consequence of the SD equation similar to the case of QED in (2+1) 
dimensions. 0* 

In our previous paper0 we stressed that the composite vector boson is actually a 
dynamical gauge boson corresponding to the hidden local symmetry (HLS) © being broken 
spontaneously, which means that the Lagrangian can be rewritten into the manifest gauge 
symmetric form. Combined with the nonlocal gauge fixing, this manifest gauge degree of 
freedom plays an important role in construction of the SD equation to keep the consistency 
between the ladder approximation and the Ward-Takahashi identity. Based on this SD 
equation, we performed the proof that the fermion dynamical mass is actually generated at 
finite value of iV (N: number of four-component fermions) in D(2 < D < 4) dimensions,0 
through the context of Ref. ^Up. Moreover we obtained explicit form of the mass function and 
a critical value of N in the limit of g — > oo in (2+1) dimensions, where g is the dimensionless 
four-fermion coupling constant of the Thirring model. 

Further studies were done by several authors. ©S'lll' The author of Ref. 13) derived the 
explicit form of the critical value of iV for finite g (denoted as N cr (g)) in the D-dimensional 
Thirring model by using the inversion method instead of the SD equation. The explicit form 

*' The original proof in Ref. |l2|) includes the case of D — 2, although the exact solution© seems to tell 
us no mass generation. In Appendix ^ we will show why this difference occurs. 
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of N cr (g) is somewhat different from the one given in Ref. [13]) and this paper (See §|4] 0), 
although they agree with each other qualitatively as to the existence of the critical N. In 
Refs. |H|),|I5D the lattice simulation were performed in (2+1) dimensions, and the results 
seem to support those of analytical studies. EH^E!) 

In this paper we first recapitulate the formalism of Ref. O) for the Thirring model as a 
gauge theory and the SD equation in the nonlocal gauge. We then give a revised proof of 
the existence of a nontrivial solution, since the previous proofS was based on an implicit 
assumption and thus was not complete. 

Next we show explicit form of the analytical solution in the limit of g — > oo in D(2 < D < 4) 
dimensions not restricted to (2+1) dimensions. By introducing a certain ansatz for the kernel 
of the SD equation, we further obtain explicit form of the analytical solution for finite value 
of g not restricted to g —>■ oo. The validity of the ansatz will be checked by the numerical 
solution in the case of (2+1) dimensions. The analytical solutions for the finite coupling 
constant will provide us with the information for investigating the phase structure, i.e., the 
critical line in (N, 1/g) plane, the beta function, the anomalous dimension, etc.. 

This paper is organized as follows. In §2 we briefly summarize the HLS and nonlocal 
gauge fixing procedure. In the first subsection of the §3 the review will be continued to 
introduce how to construct the SD equation, and in the last subsection we will improve the 
existence proof of a nontrivial solution of the SD equation as a supplement to the previous 
one.© In §4 the analytical calculations for the SD equation in D(2 < D < 4) dimensions 
will be done to obtain explicit form of the solutions used to investigate the properties of 
the model. §5 will be devoted to the comparison between the analytical and the numerical 
solutions in (2+1) dimensions. The conclusion and discussion will be given in §6. 

§2. Hidden Local Symmetry 

2.1. Hidden Local Symmetry in the Thirring Model 

In this and the next sections we would like to review briefly the formulation used in 



Ref. [12|) to investigate the dynamical mass generation of fermion in the Z)-dimensional 
Thirring model. 

Let us start with the following Lagrangian (the Thirring Model): 

N r*i N 

CtU = E - ^ (£ Val^a) 2 , (2-1) 

a=l ZiV a=l 



*' Although we give only the critical value of g, as g CI (N), we can read N cr (g) by inverting g — g cr (N cr (g)) 
with respect to N cr (g). 
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where it is supposed that each of Dirac gamma matrices is represented formally as 4 x 4 
matrix even in D(2<D<A) dimensions, and the subscript a runs over the color of fermions. 
(Hereafter we will suppress summation symbol, X) =i-) ^ * s we ^ known that this Lagrangian 
can be rewritten with the vector auxiliary field, A^, such as 

£ Thi = * a *r fe, - -hA,) V a + -Lsj, (2-2) 



^ y/N a 2G 

where it is easy to see that A^(x) is related to fermions 



K = —7=^a7^a, (2-3) 



<N 

from the equation of motion. Eq.(2~2") has no local symmetry associated with A^ due to 



the "mass term" appearing in the last of Eq.( |2-2| ). Nevertheless, we can rewrite Eq.( [2-l| ) or 
equivalently Eq.(^) so that it actually has the hidden local symmetry (HLS) 111^0 >lll> 
associated with the gauge field A^ which is introduced through the following field redefini- 
tion: 

W a {x) = u j '(x)^ a (x) = e- i ^ a (x), (2-4a) 
A^x) = iVNv^ ■D lt u = A lt -VNd^, (2-4b) 

where the field u(x) = e % ^ x > can be recognized as the nonlinearly- realized basis of the 
NG boson field <f>(x) accompanying the spontaneous symmetry breaking of the hidden local 
symmetry, £7(l)hidden- Thus, using the fields defined above, the Lagrangian ( |2-1| ) is rewritten 
into the following form, 



N ■ 

^Thi = Ipaifllpa ~ ^q( U D H U ) 



^ ~ 7n A ») 4>* + ^5 (A. - ^d,4>)\ (2-5) 



2G 

in which the HLS is realized manifestly at the Lagrangian level: 



V>« ^« = e fa Va , Ap^A^Ap + VNdpa, <p^<p' = <p + a. (2-6) 



"Gauge equivalence" between Eq.(2_2) and (2-5) can be shown by taking particular choice 



of the gauge transformation in Eq.(2_6) as a(x) = — (j)(x), called the unitary gauge. If we 
further identify ip' a = \V a and A' = A^ under such a gauge, Eq. ( |2 6| ) is identical to Eqs.(2-4). 

Why should the existence of the HLS be emphasized? The reasons are the following: 
First, coupled with the BRS formalism, this local symmetry might enable us to prove 
the S-matrix unitarity more straightforwardly. Secondly, actual calculations, particularly 
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loop calculations, are generally hopeless in the unitary gauge, while the HLS provides us to 
take the most appropriate gauge for our purpose. For these reasons, the BRS gauge fixing 
procedure should applied to Eq.(2-5) instead of the unitary gauge. Moreover, we should 
introduce a derivative- (momentum-) dependent gauge fixing parameter £(<9 2 ) (£(/c 2 )), i.e., 
the nonlocal gauge fixing^ since it is needed to guarantee the Ward-Takahashi identity in 
the ladder SD equation (see §|3p. It has also been proven that we can construct the BRS 
symmetric Lagrangian even for the nonlocal gauge. 0) 
Here we adopt the gauge fixing term such as 

£gf = -\ Ua» + ^^0) Ua" + ^^0) , (2-7) 

where £(<9 2 ) is a nonlocal (momentum dependent) parameter. Therefore the total Lagrangian 
is given by 



--total 



^Thi + ^GF — £ip,A + Cfa (2- 



£ = ^0^0-^(e(9 2 )0)0, (2-10) 



where the 4> was rescaled as y N/ G<p 1— > (p. In this gauge fixing the fictitious NG boson <p is 
completely decoupled independently whether £ is nonlocal or not. Equation (2-9) appears 
as if we added the "covariant gauge fixing term" to Eq. (|2-2j ). However, without introducing 
the HLS, there is no reason why we should add such a term. Such a confusion was actually 
made by some authors©' who happened to arrive at the Lagrangian having the same 
form as Eq. ([2-2|) in the case of constant £. Here we stress again that this Lagrangian (2-9), 
whether the gauge parameter is nonlocal or not, can only be justified through the HLS in 
the i?£ gauge. 



§3. Schwinger-Dyson Equation and Dynamical Mass Generation 

The purposes of this section are to briefly review how to construct the SD equation of the 
.D-dimensional Thirring Model and to improve the existence proof of a nontrivial solution 
of the SD equation given in Ref . [TJ|) . 
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3.1. Schwinger-Dyson Equation in the nonlocal gauge 

The SD equation for the fermion full propagator S(p) = i[A(—p 2 )p' — B(— p 2 )]" 1 , with 
B{p) being the order parameter of the chiral symmetry breaking, is written as follows: 

(A(-p 2 )-l)$-B(-p 2 ) 

1 f A(-q 2 )i + B(-q 2 ) 

= ~N J U2^^ A 2 (-q 2 )q 2 -B 2 (-q 2 ) F ^ ^ ^ " ^ ^ 

where r„(p, q) and D fMU (p — q) denote the full vertex function and the full propagator of the 
dynamical gauge boson, respectively. As indicated in our previous paper, the R.H.S. of 
Eq.(3-1) is 0(1/N), which implies that there is no nontrivial solution in the 1/N expansion. 
Thus we should solve it and find the solution which depends on 1/N in a non-analytic way 
similar to the case of the QED 3 .0) 

We should apply some appropriate approximations to Eq.(3T) so as to reduce it to the 
soluble integral equation for the mass function M(—p 2 ) = B(—p 2 )/A(—p 2 ). Approximations 
adopted here are the following: First, bare vertex (ladder) approximation is used to -T M (p, q), 
i.e., r^{p 1 q) = 7 M . Secondly, only the one-loop diagram of the bare fermion propagators iff 
is taken to evaluate the vacuum polarization tensor n^ v (—k 2 ). This approximation may be 
consistent with both A(—p 2 ) = 1 in the nonlocal gauge fixing and the bifurcation theory 
used in later sections, which reduces the SD equation as a nonlinear integral equation to the 
linear one with respect to B(—p 2 ). 

The vacuum polarization tensor in these approximations is given by 

n„„(k)=(g» v -^)n(-k% (3-2) 



U(-k 2 ) = -C-J ■ (-k 2 ) 2 ^, C-J ee j^F(2 - f )S(f , f ), (3-3) 



where tr/, equal to four by an assumption argued in §|2|, is the trace of unit matrix in spinor 
indices, and thus the gauge boson full propagator is as follows: 



iD„ v (k) = id(-k 



2\ 



k 



2 



(3-4) 



d{ k) - G ^-n{-k 2 y v{ k) ~ ^-k 2 )G-^k 2 ■ {6b) 

Then the SD equation (3-1) is reduced to the following coupled equations for A(— p 2 ) and 

B(-p 2 y. 

d D q A{-q 2 ) 



Np 2 J i{2n) D A 2 {-q 2 )q 2 - B 2 {-q 2 
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x d(-k 2 ) 
B(-p 2 ) = - 



M~k 2 ) + 2 - D}(p.q) - 2{k - p l {k - q \ (-k 2 ) 



k 2 



1 

TV 



d D q 



B{~q 2 



i(2n) D A 2 (-q 2 )q 2 - B 2 {-q 2 ) 

x d(-k 2 )[D - rj(-k 2 )], 



(3-6a) 



(3-6b) 



where k tl =p il - q^. 

It is generally difficult to deal with the coupled equations. We follow the nonlocal gauge 
proposed by Georgi et al.E^S which reduces the coupled SD equations into a single equa- 
tion for B(—p 2 ). Due to the freedom of gauge choice of £(—k 2 ), the R.H.S. of Eq.(3-6a) can 
be set to zero, which implies A(— p 2 ) = 1. In such a nonlocal gauge B(—p 2 ) itself is a mass 
function, i.e., M(—p 2 ) = £?(— p 2 ). 

Requiring A{— p 2 ) = 1 in the nonlocal gauge, Eqs.(3-6) are translated into the following 
forms (hereafter the Euclidean notation is used): 



/-7T 

Jo 



d6sin D 6 



1 



D-ldk- 



{d{k 2 ){ V {k 2 



D-2)} 

7](k 2 )d(k 2 ) 



B ^ = ^ t D 2 d(q D - 2 )K(p,q;G) 
J\ Jo 

where the kernel K(p,q; G) is given by 
K{p,q;G) : 



q 2 B{q 2 



q 2 + B 2 (q 2 )' 



(3-7a) 
(3-7b) 



2)r 



■D-l- 



d6sm D ~ 2 6 



xd(k 2 )[D - r](k 2 )}. 



p 2 + q 2 — 2pqcos9. 



(3-8a) 
(3-8b) 



The kernel Eq.(3-8a) has the following properties; symmetry under the exchange among 
first two arguments as K(p,q;G) = K(q,p;G), positivity and finiteness summarized as 
< K(p, q;G)<oo for G > 0. Here the ultraviolet momentum cut-off A is needed since 
the integral of the R.H.S. of Eq.(3-7b) diverges according to the power counting. Now the 
role of the equations (3-7) appears to be separated. Eq.(3-7a) determines the momentum 
dependence of rj{k 2 ) (£(& 2 )) as 

(D - 2) '- k2 



rj(k 2 



(3-9) 



m D - i )d{k 2 ) Jo 

which is obtained by integrating the R.H.S. with respect to k 2 . Eq.(3-7b) is therefore reduced 
to the integral equation with respect to the only one function, B{p 2 ). 
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3.2. Existence Proof of a Nontrivial Solution and the Critical Line 

From now on we reconsider the existence proof of a nontrivial solution of the SD equation 
(3 -7b), based on the bifurcation theory discussed in Refs. ^0|),^7|) and briefly summarized in 
Appendix [A]. 

Eq.(3-7b) always has a trivial solution B(p 2 ) = 0. We are often interested in the vicinity 
of the phase transition point where the nontrivial solution also starts to exist without gap. 
Such a bifurcation point is identified by the existence of an infinitesimal solution 5B(p 2 ) 
around the trivial one B{p 2 ) = 0.0* Then we obtain the linearized equation for 5B(p 2 ): 

5B(p 2 ) = i [ A ° / d( q D - 2 )K(p,q;G)5B(q 2 ), (3-10) 

where the IR cut-off m is introduced. It is enough for us to show the existence of a nontrivial 
solution of the bifurcation equation ( p-lO ).!^ Particularly, we can obtain the exact phase 
transition point where the bifurcation takes place. Since the solution is normalized as m = 
5B(m 2 ), m is nothing but the dynamically generated fermion mass. 

Rescaling p = Ax^^ and 5B(p 2 ) = AE{x), Eq. (|3-10| ) can be rewritten as follows: 

1 r 1 



Z(x) = -f dyK(x,y;g)E(y), (3-11) 

iV J Om 



where a m = (m/A) D 2 (a m < 1), the dimensionless four-fermion coupling constant g = 
GA D ~ 2 and the dimensionless kernel is defined by 

K(x,y;g) = K(x^,y^;g). (3-12) 

As we mentioned in the end of the previous subsection, the kernel K(x, y; g) is positive, 
symmetric and finite: 

< K(x,y;g) = K(y,x;g) < oo, for x, y <E [a m , 1], g > 0. (3-13) 

These are the most important properties for the existence proof of a nontrivial solution. S 
Here let us consider the linear integral equation as an eigenvalue problem: 



<p(x) = \ f dyK{x,y;g)<t>{y), (3-14) 

A J (Jm 



whose eigenvalues and eigenfunctions are denoted by X n (g, <y m ) (|A n | > |A n +i|; n = 1, 2, . . ..) 
and 4> n {x), respectively. The kernel K(x,y;g) is a symmetric one and hence satisfies the 
following property: 

J2 X l(9,o-m) = / / dxdy[K(x,y;g)} 2 < oo. (3-15) 
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The R.H.S. of Eq.( |3T5|) gives the upper bound for each eigenvalues A n (g,cr m ), so these 
are finite. Furthermore, using the positivity of the symmetric kernel (see Eq. fl3T3|) ), it is 
easily proven that the maximal eigenvalue Xi(g, a m ) is always positive and the corresponding 
eigenfunction <pi{x) has a definite sign (nodeless solution). Moreover, we can also show 
that Xi(g,a m ) is monotonically decreasing function for <r m , i.e., Xi(g,a m ) < Xi(g,a' m ) for 
o~ m > o~' m , as proven in Appendix [B| This behavior of the eigenvalue for a m was implicitly 
assumed in the previous proof. 

In the bifurcation equation (|3-11|) this implies the following: If iV is smaller than the 



maximal eigenvalue of the kernel at a m = : N < Xi(g, 0), then there exists a value of <j m 
at which N = Xi(g,a m ) and the corresponding nontrivial nodeless solution U(x) = <pi(x) 
besides a trivial one. The above statement means that the parameter a m corresponds to the 
dynamically generated mass m = A(a m )^^ . Moreover, N CT (g) introduced as 

N CT (g) = X 1 (g,a m = 0) (3-16) 

determines the critical line in (N, 1/g) plane and separating it by two parts which are cor- 
responding to the broken and the symmetric phase, since there is no non-zero solution for 
N > N a (g) by definition. Existence of the critical line, N = N cr (g) or g = g cr (N), in the 
two-parameter space is somewhat analogous to that in the gauged N JL model. & 

Finally, we can point out something more about the critical line by using the explicit 
form of the kernel. From Eq.(3-8) it can be seen that the kernel vanishes when g = 0. This 
means that Xi(g = 0, a m ) = 0, therefore the critical line starts from (0, oo). On the other 
hand, the kernel remains finite when taking g — > oo owing to the IR cut-off, which concludes 
that the critical line ends at (iV cr (oo),0) with non-zero iV cr (oo).0 



§4. Analytical Solutions for D(2< D <4) Dimensions 

Without solving the SD equation (3-7) explicitly, it has been proven that there exists the 
critical line in (N,l/g) plane, which starts from (0, oo) to (iV cr (oo), 0). However, we have 
to find an explicit form of the solution to know more about the properties of the model, 
which are, for example, momentum dependence of the mass function, the beta function 
and the anomalous dimension of the operator ipif), etc.. Therefore we attempt to solve the 
SD equation only for the particular regions of the coupling constant, since the integral for 
the kernel, Eq.(3-8a), cannot be performed for general value of g. Furthermore we restrict 

*) At this stage we cannot yet conclude that iV cr (oo) is finite since f Q dxdyK{x, y; oo) 2 = oo. However, 
the explicit calculation in the next section tells us that ./V cr (oo) is actually finite. 
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ourselves to analyze the bifurcation equation ( ft -ICQ instead of Eq.(3-7b), since it is difficult to 
solve Eq.(3-7b) being a nonlinear integral equation. Thus we will investigate the properties 
of the model near the critical line in (N, 1/g) plane, in which we are most interested. 



4.1. Analytical Solution : g — > oo case 

First we consider that g is infinity, which is identical to the vanishing limit of the dynam- 
ical gauge boson mass proportional to 1/g. In this case the spontaneously broken t/(l)hidden 
becomes a manifest symmetry of physical states, therefore we expect that the result is similar 
to that of D-dimensional QED, which we will discuss in the end of this subsection. 

The functions of the propagator of defined in Eq.(3-4) have now the following form. 

1 



d{k 2 ) 
V(k 2 ) 



-n(k 2 ) ' 

(D - 2) 2 
D ; 



C n k 2 - D , 



(4-1) 
(4-2) 



where r)(k 2 ) is determined by Eq. (|3-9|) and is independent of the momentum k 2 , which means 
that the gauge fixing parameter £(/c 2 ) is actually nonlocal according to Eq.(3-5): 



(D-l)(D-4) CDk4 . D _ 



D 



(4-3) 



Furthermore, the integral kernel of the SD equation Eq.(3-8) can be calculated exactly as 

/ 9 \°- 2 

K(p, q; oo) = K D 



+ q + \p - q\ 



K 



D 



9{p 



D-2 



y I ^D-2 



+ 



D-2 



^ ' nD-2 



K 



D 



(D-i)r(D) 



D(D-2)r(2-%)(r(f)r 



(4-4) 
(4-5) 



By using the same rescaling as Eq. (ft-ll|) , such as 



x 



(p/A) 



D-2 



(q/A) 



D-2 



Or. 



(m/A) D 2 and 2J(x) = B(p 2 )/A, we can rewrite the SD equation as 



27(x) 



dy 



(x-y)- + 9(y-x)- 
x y 



(4-6) 



As usually done for the SD equation of QED, the solution of this integral equation is obtained 
by translating Eq. (|4-6| ) into the following differential equation with the boundary conditions 
and the normalization condition: 



d_ 

dx 



x 



,dE{x) 
dx 



N°° 

4#^> 



0. 



(4-7) 
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(xS)' \ x=1 



0. 



a 



0, 

1/(13-2) 



(IR B.C.) 
(UV B.C.) 

(Normalization Condition) 



(4-8a) 
(4-8b) 
(4-8c) 



where, 



AK 



D 



A(D - i)r(D) 



(4-9) 



D(£>-2)r(2-f)(r(f)) 3 

is the critical value of A" which corresponds to AT cr (oo) appearing in the subsection 3.2., as 
will be shown later. 

The solution of Eq. 04-71) with the IR B.C. as well as the normalization condition has 
different form for each value of ./V: 

• Region I (N < JV£>): 



Six) 



Region I (N = JV; 



oo \ . 
cr / ' 



D-2 
(J m 



sin (|(5) 



a 



in 



2 [ UJ 

sin < — 

x J 12 



i X r 

In h 6 



'V°° 

— -1, 5 = 2w- 1 tan- 1 cj. 
V 



(4-10) 



Region H (JV > A^ 



to 





1, 2 

- In hi 




-2 cr m 



1 

-D-2 



— 1 2 sinh i — 
sinh U/ Sm 1 2 



In — + 6' 

f 771 



TV ' 



5' = 2cy- 1 tanh-V. 



(4-11) 



(4-12) 



Since only the solution of Region I can satisfy the UV B.C. for a m ^ 0, we conclude that 
the dynamical mass is generated only in case of N < N™, which means is actually the 
critical value separating the symmetric phase and the broken one. As already mentioned in 
our previous paper, in (2+1) dimensions the critical value = 128/37T 2 is identical to 
the one in QED3.EI) Furthermore, the UV B.C. Eq.(4-8b) gives the relation between A" and 
'm/AY D -V as 

1,1 r 1 i 

(4-13) 



10 
2 



In — + 25 



1,2,-.. , 
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where the solution with n = 1 corresponds to the nodeless (ground state) solution whose 
scaling behavior is read from Eq.( |4-13| ): 



m 
A 



26 
6 D~2 



exp 



2tt 



{D-2)JN~/N-1 



(4-14) 



As was done in Refs. ^),|T2]), we recognize it as determining the cut-off dependence of 
iV = N(A), by which we can have a finite m when A is taken to be infinite (continuum 
limit). In fact, if N(A) goes to continuously when A — > oo, it is possible to keep m 
finite, and hence the beta function which determines the behavior of N(A) with respect to 
A is derived as 



Pa(N) = A 



dN(A) (D 



-N 2 



cr 

N 



1 



(4-15) 



dA irN™ 

where this form is valid for N being nearly below N™. Eq. Q4-15 ) implies that the critical value 
is a UV fixed point on which the continuum theory will be determined. Furthermore the 
anomalous dimension of the operator i/ji/j, which is relatively large (7^ ~ 0(1)) in models 
inducing the dynamical mass generation, is obtained by 

d\n{$ip) A 



7^ (AO 



din A 
D-2 

— 4 

D-2 



2N(3 A {N) 



1 - 



AN 3 /N: 



cr 



N 



- 1 



(4-16) 



In Eq.(4-16) the vacuum expectation value (ipip)^ is given as 

32 



• A^ 2 • m 2 A 2 



(47r)T(D-2)r(f)JV2 
where we used the following expression for (ipip)A in the bifurcation theory: 



(4-17) 



2tr/ N 



(4vr)fr(f; 



dp p 



D-l 



B{p) 



p 2 + B 2 (p) 



2tr/ N 



D N 



dpp D - 3 B(p), 



(4-18) 



and substituted B(p) into the solution in Region I. As is expected, Eq.(4-16) tells us that 
the anomalous dimension 7^(iV) is actually large, especially ^^(N) ~ (D — 2)/2 in the 
vicinity of the critical point. 

The above renormalization group functions are obtained along with the horizontal line 
in (N, 1/g) plane, which is specified by 1/g = 0. However, it might not be natural to give 
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the cut-off dependence to N which should be an integer number. Thus we hope to know the 
solution for a finite g to attach the yl-dependence not to N but to g, although we cannot 
solve the SD equation exactly for any g. In the rest of this section we will solve the SD 
equation Eq.( ^-10[ ) for finite g though restricted region. 



4.2. Analytical Solution : g 3> 1 case 

First we attempt to find the next-to-leading order solution in 1/g expansion, which means 
that we expands Eq. ( |3- 10|) in 1/g up to 0(1/ g). In this expansion the integrand of Eq.(3-8) 
is reduced to 

(D-2f 



V (k 2 



D 



{l - (D - ^CdG- 1 ^- ) + 0(G~ 



2un /,2M KD-1)C D 1 



d(k )[D - g(k )] 



where 



D 



D k D ~ 2 + f D /G 

_ (4 -£>)£> 



0(G~ 



-C 



D 



(4-19) 
(4-20) 

(4-21) 



is a value of 0(1) in 2<_D<4. Eq.(4-20) seems to have a simple form, but we use the ansatz 

k 2 = p 2 + q 2 — 2pq cos 6 ~ max (p 2 , q 2 ^j (4-22) 

for Eq.(3-8b) to calculate the integral kernel analytically. This ansatz will be checked by the 
numerical study in §[| even in (2+1) dimensions. By substituting Eq.( |4-22| ) into Eq.(4-20), 
we can evaluate Eq.(3-8a) as 

1 



K(p,q;G) = K D 



6(p 



D-2 



Q 



D-2\ 



+ (p <-> q) 



(4-23) 



J p D ~ 2 + fo/G 

where Kp is defined by Eq.(4-5). Similarly to the previous subsection, by rescaling with A 
and translating into the differential equation, we obtain 

d ( , , . . 2 dE(x) 
cTA {X + fD/g) —l ' AN 



N°° 



0. 



(4-24) 



E'(a m ) = 0, (IR B.C.) 

[(x + f D /g)S]' x=1 = 0, (UVB.C0 



(4-25a) 
(4-25b) 



with the normalization condition being the same as Eq.(4-8c). The solution of Eq. (|4-24|) 
is easily obtained by substituting x into x + fn/g in the solution of g — » oo case, and the 
explicit form of the solution is given by 



E(x) 



~ D-2 



'o- m + f D /g^ 



sin (f 5) \ x + Id/ 9 



sin 



In X + fD J 9 +5 



0~ n 



Jul 9 



(4-26) 
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In Eq. ( |4-26| ) to and S are the same as those defined in Eq. ( f4-10|) , and we wrote down only the 
solution corresponding to the case of N < N™, because only this solution has a possibility 
to satisfy the UV boundary condition Eq.(4-25b). The UV boundary condition tells us the 
scaling relation 



"' {(l + f»/9)e 2S -%-fn/g} D - 2 =(l + gJfn)^ 



1 _ h 



(4-27) 



A 

where, 

g cr = f D (e^ 25 - l) (4-28) 



is the critical coupling above which the dynamical mass is generated. From Eq.( 4-27j) the 
beta function of g can be seen as follows: 

Mg ) = A a JL = { D- 2) g (l - JL) , (4-29) 

which coincides with the one of the D-dimensional GN model!)' B except for the expression 
of g cr . At this moment it is not clear why these are the same to each other in contrast to the 
case of g <C 1 to be discussed in the next subsection. Furthermore the anomalous dimension 
is given as 

1 ^{g) = (D-2) \ + J D,9 ;\ (4-30) 
1 + 3D 1 9 

which is calculated by using 

39/V 2 

(M A = 5 (a m + f D /g)mA D ~ 2 (4-31) 

obtained by substituting Eq. (|4-26Q into Eq.(4-18). Again 7^,(5 ^ g cr ) — D — 2 in coincidence 
with GN model in D dimensions. 

4.3. Analytical solution : g -C 1 case 

We finally analyze the case that the four fermion coupling constant is very small. In this 
case the integral defining the kernel function of the SD equation can easily be performed 
explicitly as 

d(k 2 ) [D - r](k 2 )} =G + 0(G 2 ) } (4-32) 

2DG 

K( *< G > = m^HD - mm s MdG - (433) 
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which reduces a momentum independent soluition. Then the scaling relation is obtained as 

l/(D-2) 

where 



m ' 1 - — I , (4-34) 



iV (47r) g / 2 (D - 2)r(D/2)N 

= Md = 2D (4 ' 35) 

is the critical coupling in this case, and the beta function and the anomalous dimension are 
obtained as 

A (g) = (D- 2)g f 1 - — ) , (4-36) 
j^(g) = (D-2)-^, (4-37) 

9cT 

respectively, similarly to that of the D-dimensional GN model. i^'B 1 

As a summary of this section let us notice that both g ^> 1 and g <C 1 cases give the 
same result 7^ ~ (D — 2) at g ~ # cr , which is expected from the fact that the anomalous 
dimension near the critical line is a kind of the critical exponents. On the other hand, the 
result obtained in the infinite limit of the coupling constant is different from the others, 
since in this case the renormalization group functions are calculated along with the A^-axis 
in (AT, 1/g) plane. The similar situation can be seen in the gauged NJL model. 0> 



§5. Numerical Study for the SD equation in (2+1) Dimensions 

As indicated in the last section, we check the validity of the ansatz used to solve the SD 
equation for g ^> 1. 

In this paper we would like to show the result only in (2+1) dimensions, in which the 
integral kernel can be calculated for any value of g as 



where, 



K(x,y;g) 



K{z;g) = - 
9 



TT 2 xy 



(x +y,g)- k(\x - y\;g)] 



Agz 1 / 1 

3 gz V g 2 z\ 



ln(l +gz) 



(5-1) 



(5-2) 



In the above equations, we already used variables normalized by the cut-off, which are used 
frequently in previous sections. Since the kernel is solved analytically, we can skip a few steps 
in the numerical calculation, which include complicated numerical integrations. Moreover, 
let us restrict to the linearized SD equation Eq. (|3-10|) instead of the original nonlinear integral 
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equation, since the linear integral equation can be treated as an eigenvalue problem which is 
easier to treat on computational analysis. As mentioned in §^|, it is enough for investigating 
the structure of the critical line. 

The figures following below show the result of the numerical calculation with the analyt- 
ical solutions given in the previous section. 




Fig. 1. The critical line in (2+1) dimensions 

The prot points in Fig.[l| are the result of numerical calculation with the kernel of Eq. ( |5T|) . 
Lines represent the analytical solutions, a solid one is the solution of g 3> 1 given in §4.2., 
and a dashed one is that of g <C 1 given in §4.3.. 

From these figures it is concluded that that the analytical calculation with the ansatz 
well reproduces the critical line obtained by the numerical analysis. 

§6. Conclusion and Discussion 

As an extension of our previous work,© we have further studied dynamical fermion 
mass generation of the Thirring model as a gauge theory in D(2<D <4) dimensions. 

After the introduction of the hidden local symmetry and the Schwinger-Dyson equation, 
we have completed the existence proof of a nontrivial solution of the SD equation made in 
Ref. |12|) where an assumption on the property of eigenvalue of the integral equation was 
implicitly made. Furthermore we have shown that the critical line in (N, 1/g) plane starts 
from (0, oo) and ends at (iV cr (oo), 0), which can be proven only by using the asymptotic 
behavior of the kernel with respect to g, the dimensionless four-fermion coupling constant. 

As natural but nontrivial extensions of the previous work, we have found the analytical 
solutions of the linearized SD equation (bifurcation theory) in D(2 < D <4) dimensions for 
some region of g. 
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In the case of g — > oo, the SD equation in D dimensions can be solved exactly and the 
nontrivial solution have been obtained in the region N < N^(< oo), which is consistent 
with the result in (2+1) dimensions. S 1 This general D-dimensional form of the solution 
tells us that the critical value which depends only on D is monotonically decreasing for 
2 < D < 4 from infinity (D — > 2) to zero (D — > 4), and reproduces the result in D = 3 of 
Ref. O). This result would imply that in D = 4 dynamical mass generation does not occur 
even if we takes g — > oo, while there is only the broken phase in D = 2. However, there are 
subtleties in both dimensions concerning the UV/IR divergences which we did not take into 
account because of the absence of those in 2 < D < 4 dimensions. As indicated in Appendix 
[A], we should not believe the result of the linearized SD equation in D = 2, although the 
SD equation itself is exact for D = 2, since the IR singularity may break the validity of 
this equation. In case of D = 4 the UV divergence should be properly regularized in the 
evaluation of the propagator of the dynamical gauge boson. Therefore we should study more 
carefully in this case. 

The analytical solutions have been obtained also for finite g in the cases of g ^> 1 and 
g C 1. These have given explicit form of the critical line which indicates a finite critical 



coupling constant g cr for iV < as opposed to the result of Ref. 11) claims g CI for all 



N. On the other hand, at least qualitatively, our result seems to be consistent with that of 



Ref. [13]) which is obtained also based on the HLS but through the inversion method instead 
of the SD equation. Furthermore, it should be emphasized that the numerical studies0'@ 
would indicate that there is a finite critical value of N. 

Moreover, we have found that the solutions obtained at both g ^> 1 and g<l gives the 
same expression for /3A{g) and the same value of the anomalous dimension on the critical 
line, i.e., 7 m (g) = (D — 2) at g — g CT which is identical to that of D-dimensional Gross-Neveu 
model. & 

Finally, we have confirmed validity of the ansatz used in the analytical calculation through 
numerical study of the linearized SD equation with the exact kernel, although limited to 
(2+1) dimensions. This result has shown that the numerical solution agrees with the analyt- 
ical one in both qualitative and quantitative sense. However, it may be nontrivial whether 
or not we can extend this agreement to other dimensions, which should be confirmed by the 
efforts to the more complicated numerical studies. 
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Appendix A 

Brief Summary of The Bifurcation Theory 



In this Appendix we will briefly summarize the bifurcation theory, and comment on a 
possibility of break down of the bifurcation theory in (1+1) dimensions. 
The bifurcation theory is applied to such a following functional equation, 

4>(x) = XF[4 > ;x), (A-l) 

where, F[<f>] x) in the R.H.S. is a functional of 4>(x) as well as a function of x, and A is a 
parameter. Here let us assume that one of the solutions of Eq. ( |AT ) is already known, which 



is denoted as <f>o(x; A). Hence we are interested in the solution which bifurcates from 4>o(x; A) 
continuously at A = Ab r - 

As discussed in Ref. ^9|), the necessary condition for the existence of such a bifurcated 
solution is as follows. The functional F[4>; x) is functionally differentiate in the third order 
with respect to (f>(x) at <p(x) = 4>o(x; Abr) , where the bifurcation point Ab r is determined as 
the eigenvalue of the linearized equation for 5<p(x): 



5<f ) (x) = Xld y 6mx) 



my) 



8(f>(y). (A-2) 

{x)=<f> (x) 



Thus the bifurcation solution is given as <f>(x) = 4>o(x; Ab r ) + 6(f)(x) in the vicinity of Ab r - 

We can easily check with the SD equation (3-7b) that the R.H.S. satisfies the necessary 
condition. Moreover, we can also confirm that the necessary condition is satisfied even if we 
use the fermion full propagator when calculating the dynamical gauge boson's propagator 
and the integral kernel of Eq.(3-7b), although the kernel cannot be solved explicitly in such 
a case. 

Concerning with the above statement, we should comment on the case of D = 2 in 
which the exact solution was found to be equivalent to a free theory having no dimensionful 
parameter. However, if we solve Eq.( p-10|) in D = 2 naively, a nontrivial solution for g > 
exists, which is constant for the momentum as 5B(p 2 ) = m — Aexp{— N(l + n/g)}. It 
should be emphasized that the ladder approximation is not wrong in this case, although it 
has been usually suspected. The reason is that the set of vertex functions 

S(p)=i/h D ia ,(k) = (g- 1 + Tr- 1 )g ia , i r>(p,q)=>y», 
other vertex function = 

is the exact solution of the full series of the SD equations, which means that the ladder 
approximation is exact in (1+1) dimensions, if the bifurcation theory is valid. 
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Do these considerations conclude the dynamical mass generation in the (l+l)-dimensional 
Thirring model? This confliction will be resolved by investigating the integral kernel more 
in detail. Thus it might be shown that we are not able to apply the bifurcation theory to the 
case of D = 2 because of the non-analyticity of K(p, q; G) with respect to the mass function 
B{p 2 ) near a trivial solution B{j> 2 ) = 0, i.e., 5K(p, q; G) / 5B(j> 2 )\ b ^=q = oo. This property 
of the kernel is contrary to the necessary condition explained above. Therefore in (1+1) 
dimension we should not naively conclude the existence of a nontrivial solution through the 
bifurcation theory. 



Appendix B 

Proof that X(g,a m )is a decreasing function of o~ r , 



In §^ we needed to prove that the maximal eigenvalue of the integral equation ( 3-14]) is 



monotonically decreasing function with respect to a m . For this purpose, we will prove more 
general statement in this Appendix. 

Let us consider the following eigenvalue problem: 

<Kx) = \ ! dyK(x,y)<Ky), (B-l) 



A J a 

where K(x, y) is a real symmetric integral kernel which is finite on the region a < x,y < b. 
We denote eigenvalues and corresponding eigenfunctions as A n (a) and 4> n (x;a), in which n 



is assigned as Ai(a) > A2(a) > .... The dependence of a is defined by both Eq. (|B-l|) and the 
normalization condition 

b 

dx(f) m (x] a)(f) n (x] a) = 5 mn , (B-2) 



b r b 



and hence A n (a) is expressed as 

K{a)=i I dxdy(j) n (x;a)K(x,y)(f) n (y;a). (B-3) 

J a J a 

By using this explicit expression of A„(a), we can calculate the derivative of X n (a) with 
respect to a as 

dXJa) r b 

— = -2(p n (a;a) dyK(a,y)(p n {y; a) 

CLQj J a 

+ 2 f f dxdy 9< ^ n X,a K(x,y)<f> n (y;a) 

J a J a UCL 

s2 , / \ f b i d(j) n (x;a) 



-2A n (a) (0„(a; a)) +2A n (a)/ dx — (j> n (x;a), (B-4) 
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where, we use the symmetry property of K(x, y). The second term of the last line of Eq.(B-4) 
is derived by the differentiation of Eg. ( |B-2| ): 

= — / dxcf) n (x; a)<j) n (x; a) 



da J a 

= - (0 n (a; a)) 2 + 2 / dx ^jf' n (x; a). (B-5) 
Ja aa 

Substituting Eq.(B-5) into Eq.(B-4), we finally obtain the formula as 

d\n{a) w w , , ,v 2 / < (A n (a) > 0) 

— = -Anfl Woo I B-6 

aa 1 > (^n(a) < 0) 

which is the statement that should be proven. 

Combining with the fact that the maximal eigenvalue Ai(a) has always positive sign, we 
have actually proven that Ai(a) is monotonically decreasing function for the lower bound of 
the integral. 
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